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A star orbiting a supermassive black hole can be tidally disrupted if the black hole's gravitational 
tidal field exceeds the star's self gravity at pericenter. Some of this stellar tidal debris can become 
gravitationally bound to the black hole, leading to a bright electromagnetic flare with bolometric 
luminosity proportional to the rate at which material falls back to pericenter. In the Newtonian 
limit, this flare will have a light curve that scales as t~^^^ if the tidal debris has a flat distribution 
in binding energy. We investigate the time dependence of the black-hole mass accretion rate when 
tidal disruption occurs close enough the black hole that relativistic effects are significant. We find 
that for orbits with pericenters comparable to the radius of the marginally bound circular orbit, 
relativistic effects can double the peak accretion rate and halve the time it takes to reach this peak 
accretion rate. The accretion rate depends on both the magnitude of the black-hole spin and its 
orientation with respect to the stellar orbit; for orbits with a given pericenter radius in Boyer- 
Lindquist coordinates, a maximal black- hole spin anti-aligned with the orbital angular momentum 
leads to the largest peak accretion rate. 



I. INTRODUCTION 

Active galactic nuclei (AGN) are believed to be pow- 
ered by accretion onto compact objects with masses 
M > 10^ Mq [T]; such objects will inevitably collapse 
into supermassive black holes (SBHs) on short timescales 
[2]. Early work OS] conjectured that these AGN were 
fueled by the tidal disruption of stars passing too close 
to the SBHs, and that SBHs could grow to their ob- 
served masses by accreting debris from such tidal disrup- 
tion events (TDEs). Although SBHs are now believed 
to grow primarily by accreting gas driven into galactic 
centers by tidal torques during galactic mergers [5"^ 'B^, 
interest in TDEs was renewed when it was realized that 
they could power bright flares lasting several years in 
otherwise quiescent galaxies f7^. 

The Roentgensatellit (ROSAT) observed such a flare 
in soft X-rays in NGC 5905, a galaxy with no previous 
indication of nuclear activity [8 . A systematic survey 
for X-ray flares within the ROSAT All-Sky Survey dis- 
covered five such events, implying a rate of 9.1 x 10~^ 
galaxy"^ yT~^ consistent with the predicted rate of TDEs 
[9]. Five additional TDE candidates were discovered in 
the X-ray by the XMM-Newton Slew Survey [10^. TDEs 
emit in the UV and optical as well; several TDE candi- 
dates discovered in the UV by the Galaxy Evolution Ex- 
plorer (GALEX) were found to have optical counterparts 
in the Canada-France-Hawaii Telescope Legacy Survey 
and Panoramic Survey Telescope and Rapid Response 
System (Pan-STARRS) Medium Deep Survey fTTHIS]. 
Two new TDE candidates were identified through a se- 
ries of rigorous cuts on the many optical transients found 
in Stripe 82 of the Sloan Digital Sky Survey (SDSS); this 
implies that deeper, higher-cadence surveys such as that 
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which will be undertaken by the Large Synoptic Survey 
Telescope (LSST) could find thousands of TDEs per year 
[M]. TDEs can also be observed in hard X-rays if they 
launch jets pointed towards us; the events Sw J1644+57 
[15-17 and Sw J2058+05 [18 discovered by the Burst 
Alert Telescope on the Swift satellite are conjectured to 
be such relativistic tidal disruption flares. 

One of the most exciting applications of observing 
TDEs is the possibility of measuring SBH spins. The 
magnitude of a SBH's spin depends sensitively on the 
manner in which that SBH was assembled. A nonspin- 
ning SBH can attain a dimensionless spin a/M 0.998 
after increasing its mass by a factor of \/6 through the 
coherent accretion of gas in the equatorial plane [191 120] • 
If a SBH grows instead through the chaotic accretion 
of misaligned clumps of gas [2TH23] , its spin magnitude 
could be much more modest {a/M < 0.2). Studies of 
SBH spins in the context of hierarchical galaxy formation 
indicate that a wide range of distributions are possible, 
depending on one's assumptions about the relative con- 
tributions of coherent accretion, chaotic accretion, and 
binary mergers to SBH spin evolution [24-26 . Simula- 
tions in which gas is accreted coherently or chaotically 
lead to dimensionless spin distributions sharply peaked 
about unity or zero respectively at all redshifts. Real 
SBHs probably grow through an as yet undetermined 
combination of coherent and chaotic accretion, suggest- 
ing that almost any distribution of SBH spins is possible. 
This considerable theoretical uncertainty emphasizes the 
urgent need for greater observational constraints. 

Fortunately, there has recently been tremendous 
progress measuring SBH spins using X-ray spectroscopy. 
If an AGN accretion disk is surrounded by a hot, nonther- 
mal corona, a portion of the hard X-rays emitted by this 
corona may be reprocessed by the colder disk [27l [28] , 
This hard X-ray irradiation will cause the 6.4 keV iron 
Ka line in the disk to fluoresce, with much of the flux 
in this line coming from the innermost regions around 
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the SBH [291 EO]. The shape and variabihty of the hne 
profile is sensitive to general-relativistic effects including 
SBH spin. Observations of fluorescent iron Ka lines with 
the XMM- Newton and Suzaku satellites have been used 
to constrain the spin magnitudes of eight nearby AGN, 
three of which are near maximal (a/M > 0.98) while the 
remaining five are intermediate {a/M 0.7) [31] |32]. 

Despite these impressive new results, this approach to 
measuring SBH spins is inherently limited to highly ac- 
creting AGN whose spins may not reflect the SBH popu- 
lation as a whole. TDEs provide an opportunity to mea- 
sure the spins of a less biased SBH population, and spin 
estimates based on TDEs should have different system- 
atic errors than previous approaches. If a TDE launches 
jets that extract rotational energy from the SBH via the 
Blandford-Znajek mechanism [33 , the observed peak X- 
ray luminosity of the TDE can be used to constrain the 
SBH spin magnitude [34 . If the jet is misaligned with 
the SBH spin, the Lense-Thirring effect [35 will cause the 
jet axis to precess, modulating the observed X-ray emis- 
sion [36l [37] . Both of these approaches have been applied 
to the two Swift events J1644+57 and J2058+05. Mean 
SBH spins may also be estimated from the observed TDE 
rate [SSHlO]. 

In this paper, we focus on yet another way in which 
TDEs can be used to constrain SBH spins, through the 
observed light curve of individual events. Several previ- 
ous authors have investigated relativistic effects during 
tidal disruption [4TH44) : Laguna et al [42 provide light 
curves for TDEs by non-spinning SBHs. If a star with 
mass and radius is on an initially parabolic orbit 
with pericenter r^, after tidal disruption roughly half of 
the stellar mass will become gravitationally bound to the 
SBH with specific binding energy |0 
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where /3 = rud/rp is the penetration factor, rtid = 
{M /m^Y^^R^ is the tidal radius, and the binding energy 
is defined as the positive energy required to unbind the 
system. Newtonian orbits with specific binding energy E 
have orbital period 
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implying that the debris will return to pericenter a time 
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after tidal disruption [45]. Internal shocks during this 
pericenter passage will cause the debris to settle into an 
axisymmetric torus on this timescale. If the viscosity is 
high enough to allow accretion on a similar or shorter 
timescale, Eq. ([2| implies that the TDE luminosity L oc 



\dE/dt\ oc t~^^^ provided the energy distribution dm/dE 
of the tidal debris is constant ^ • The light curves of 
several TDE candidates have been reasonably well fit by 
^-5/3 pQ^gj. lawg the optical, UV, and X-ray [121 [131 

[ismi]. 

The L oc t-^/^ dependence of TDE light curves rests 
on two assumptions: the constancy of dm/dE and the va- 
lidity of the Newtonian relation ([2| between the specific 
binding energy and orbital period. The first of these as- 
sumptions was investigated by Lodato, King, and Pringle 
[48 (hereafter LKP09), who found that dm/dE did in- 
deed depend on the adiabatic index 7 of the tidally dis- 
rupted star. This dependence on 7 is reflected in the ear- 
liest portions of the TDE light curve; Gezari et al. [13] 
showed that the TDE candidate PSl-lOjh discovered by 
Pan- STARRS and GALEX is more consistent with a fully 
convective star or degenerate core with 7 = 5/3 than a 
solar- type star with 7 = 4/3. 

The second assumption was investigated by Haas et 
al. [49 , who performed six numerical-relativity simula- 
tions of the tidal disruption of a white dwarf by a 10^ Mq 
intermediate-mass black hole (IMBH). The IMBH was 
nonspinning in the first simulation, and in the remain- 
ing five simulations had a spin with the same magnitude 
(a/M = 0.6) but different orientations with respect to 
the orbital angular momentum and initial position of the 
white dwarf. They found that the prompt accretion rate 
at t < 6 s after tidal disruption could vary by almost 
two orders of magnitude as a function of spin direction. 
Numerical-relativity simulations are computationally ex- 
pensive however, and it is not feasible to perform such 
simulations for the long timescales ttid given by Eq. (|3| 
over which debris continues to fall back to pericenter. 
Haas et al. [49] instead estimate the fallback time for 
each fluid element at the end of their simulations as- 
suming that subsequent motion is on Keplerian ballis- 
tic trajectories [50]. They find, similar to LKP09, that 
dm/dt oc at sufficiently late times. 

In this paper, we use the assumption that tidal debris 
moves on Kerr geodesies to calculate mass accretion rates 
as a function of time throughout the full fallback regime. 
Our focus will be on the most massive SBHs where the 
tidal radius rtid is comparable to the Schwarzschild ra- 
dius Ts] the extreme mass ratios {m^/M < 10~^) imply 
that full numerical relativity is not required leading to 
a vast reduction in computational expense. In Sec. |ll| 
we review the Newtonian calculation of the mass accre- 
tion rate given by LKP09. We then propose a relativistic 
generalization of this calculation in Sec. HI Using this 



new relativistic framework, we present the accretion rate 
as a function of SBH mass, spin, and initial stellar or- 
bit in Sec. 



IV 



We discuss the implications of our results 

and future applications of this framework in Sec. [Vj For 
the reader's convenience, we have summarized the Kerr 
metric and the use of Fermi normal coordinates in Ap- 
pendices [A] and [Bj 
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II. NEWTONIAN ACCRETION RATE 

A simplified model of stellar tidal disruption in the 
Newtonian limit is presented in Sec. 2 of LKP09. This 
model is premised on the existence of a clear hierarchy 
between the specific orbital kinetic energy E^rh of the 
initial star, the energy £^tid gained or lost during tidal 
disruption, and the self-binding energy £^bin of the star 
[50] , These energy scales are given by 

Eovh = , 4a) 

^tid = = /3g'/'^orb , (4b) 

Eun = ^=r^g^/Xrb, (4C) 
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where Tp is the pericenter and q = m^/M <C 1 is the 
mass ratio. This hierarchy implies that to zeroth order 
in g, the specific binding energy of a fluid element of the 
tidal debris is the same as that of the center of mass of 
the initial star. At 0(g-^/^), a fluid element located at r 
receives a correction to its specific energy 
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where points from the black hole to the center of the 
star at pericenter, and vectors with hats above them are 
unit vectors. This correction is of order £^tid- The kinetic 
energy of the fluid element in the star's center-of-mass 
frame is of order £^bin, as is the gravitational binding en- 
ergy of the fluid element to the star. These can therefore 
be safely neglected to this order in q. 

Given this hierarchy of energy scales, Rees [7 approx- 
imated the distribution in specific energy of the tidal de- 
bris of a star on an initially parabolic orbit {E — 0) as 
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As described in Secjlj this approximation leads to a mass 
accretion rate that begins a time ttid after tidal disruption 
given by Eq. (|3| and falls of as t~^l'^ thereafter. 

LKP09 move beyond this level of approximation, ex- 
plicitly calculating the energy distribution of the tidal de- 
bris to 0(g^/^) for different polytropic equations of state 
P = Cp^ . They begin by solving the Lane-Emden equa- 
tion to determine the density profile of the initial star 

m 



= -9" 



(7) 



where n = (7 — 1) ^ is the polytropic index and O and ^ 
are dimensionless variables which give the star's density 



FIG. 1: The distribution dm/dE of the specific binding energy 
E for the tidal debris of a star on an initially parabolic orbit. 
The solid black curve corresponds to a adiabatic index 7 = 
5/3 appropriate for a fully convective star, while the dashed 
blue and dotted red curves correspond to 7 = 1.4 and 4/3 
respectively. 



profile through the relations 

P = PcO" , 



1/2 



(8a) 
(8b) 

(8c) 



where pc is the star's central density. The Lane-Emden 
equation ^ can be solved numerically by beginning at 
the star's center ^ = where 6(0) = 1, dG/d^{0) = 0, 
then integrating to the star's surface at ^1 defined by the 
condition 9(^i) = 0. If the specific binding energy of 
a fluid element is related to its position by Eq. ([5|, the 
energy distribution dm/dE of the tidal debris is given in 
terms of the star's initial density profile 0(0 by 



£;tid dm _ lE/Etid Q""^^^ 



m* dE 2 /q 6^x2 
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where x = ^/^i is a rescaled radial coordinate. This di- 
mensionless energy distribution for several adiabatic in- 
dices 7 is shown in Fig. [l] which reproduces the left panel 
of Fig. 2 of LKP09. 

If the orbital period t of a fluid element is given in 
terms of its specific binding energy by Eq. ([2|, the rate 
dm/dt at which mass falls back to pericenter and is sub- 
sequently accreted is given by 



to dm 
m^ dt 



1 Vr, 



Etid dm 
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where = 27r(r^/GM)^/2 is the period of an orbit with 
semi- major axis Tp. We plot this dimensionless accretion 
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FIG. 2: The rate dm/dt at which the debris of a tidally dis- 
rupted star falls back to pericenter and is subsequently ac- 
creted. The time t is given in units of to = 27r(rp/GM)^/^, 
the period of an orbit with semi-major axis Tp. As in Fig. [l] 
the solid black, dashed blue, and dotted red curves correspond 
to adiabatic indices 7 = 5/3, 1.4, and 4/3 respectively. The 
dot-dashed green line shows the canonical t~^^^ dependence 
predicted bv (71146). 



bound orbits about a Newtonian point mass. Since we 
are interested in the rate at which tidal debris returns to 
pericenter, we will focus on the radial period and ignore 
that the second pericenter will generally occur at differ- 
ent values of and ^ than the pericenter at which the 
star was initially disrupted. This precession may affect 
the timescale on which the tidal debris settles into an ac- 
cretion disk, but this process is beyond the scope of the 
current work. 

A second difference between Kerr geodesies and New- 
tonian orbits is that the radial period depends not only 
on the specific energy as in Eq. but also on the 
specific angular momentum and Carter constant Q 
as can be seen from Eq. ( |AlQb ). Determining the ac- 
cretion rate dm/dt will therefore require not just the en- 
ergy distribution dm/dE, but the full joint distribution 
d^m/ dEdLzdQ for all three constants of motion. Also, 
as we do not have a nice relation like Eq. ( 2| fo r the ra- 
dial period, we will have to integrate Eqs. (AlO) directly 
for each fluid element of the tidal debris. 

The greatest difference between our relativistic model 
and the LKP09 model is the manner in which we deter- 
mine the distribution d^m/dEdLzdQ. In the Newtonian 
limit, the specific energy 
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rate for = 100 and several values of 7 in Fig. [2j a 
reproduction of the middle panel of Fig. 2 of LKP09. 

After introducing this simplified model of tidal disrup- 
tion in Sec. 2 of their paper, LKP09 devote the rest 
of the paper to testing its validity using a series of 
non-relativistic smoothed-particle- hydrodynamics (SPH) 
simulations. They find that Eq. ([9| accurately describes 
the energy distribution of the tidal debris after the ini- 
tial density profile of the star has been homologously ex- 
panded by a factor of 1.6 to 2.5 depending on the adi- 
abatic index 7. This expansion reflects the stretching 
of the star by tidal forces prior to its arrival at peri- 
center. The mass accretion rate dm/dt measured in the 
simulations also confirms the predictions shown in Fig. |2j 
a stiffer equation of state (larger 7) leads to a higher 
peak accretion rate and a faster approach to the late- 
time behavior. The success of the LKP09 model 
provides a strong motivation for developing a relativis- 
tic version of this model appropriate to TDEs where the 
pericenter distance is comparable to the Schwarzschild 
radius rs = 2GM/c^ of the SBH. 



III. RELATIVISTIC ACCRETION RATE 

Kerr geodesies differ from Newtonian orbits in several 
important respects that need to be incorporated into our 
new relativistic model for tidal disruption. We provide a 
brief review of Kerr geodesies in Appendix |Aj One differ- 
ence is that Kerr geodesies generally do not close, unlike 



separates into three distinct components: the rest-mass 
energy c^, the kinetic energy i'^/2, and the gravitational- 
potential energy <l>(r). The rest-mass energy is neglected 
in the LKP09 model because it does not affect the orbital 
motion in the Newtonian limit. Although the kinetic 
energy of the center of mass is equal in magnitude to the 
gravitational-potential energy for a parabolic orbit, the 
relative velocity v between fluid elements at pericenter 
can be neglected because it leads to energy corrections of 
order £^bin, which according to Eq. Q are suppressed by 
a factor q^l'^ / 0^ compared to the dominant corrections of 
order £^tid arising from the gravitational potential. The 
constants of motion E^ L^, and Q given by Eqs. ( A5 ) and 



( A8 ) do not neatly separate into components that depend 



only on r or as in Eq. ( 11 ), however the recognition that 
we can set the relative velocity v = provides the key for 
calculating the relativistic distribution d^m/ dEdLzdQ. 

It is often convenient to use Fermi normal coordinates 
when calculating quantities near a point (such as the cen- 
ter of mass of a tidally disrupting star) that moves along 
a Kerr geodesic. We review the use of these coordinates 
in Appendix |B] If a fluid element is located at a fixed 
position X* in Fermi normal coordinates, its relative ve- 
locity with respect to the origin (the center of mass) will 
vanish [52^: 



A^V;3(X'Af ) = XU^V^jAf = 



(12) 



The final equality follows because our basis vectors 
were constructed to be parallel transported along the cen- 
tral geodesic (Aq V^Af = 0). The difference in specific 
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energy between this fluid element and one located at the 
origin is given by 



7^0 



dt 



(13) 



The secon d equality follows from the deflnition of E in 
Eq. (A5a). The third equality follows from the com- 



patibility of the metric with the derivative operator, 
the commutation of coordinate vector flelds (Af Vq,Aq = 
V^Af ) [52 , and the parallel transport of the basis vec- 
tors. The flnal equality follows from the deflnition of the 
Christoffel symbols. Similar calculations show that 



AK = X'XfVc^K 



a-pf 



and 



rAIl 



(14) 



(15) 



where K is deflned in Eq. (A8) and P and are given 



by Eq. (A7). We can then compute the difference 
AQ = AK- 2{L^ - aE){AL^ - aAE) 



(16) 



in the Carter constant Q of the fluid element compared 
to the center of mass. 

Our prescription for calculating the mass fallback rate 
in the relativist ic limit is as follows: 

(1) Choose an orbit for the tidally disrupting star char- 
acterized by its 4- velocity Aq at pericenter and con- 
stants of motion E^ Lz^ and Q. 

(2) Determine the density proflle p{r) of the star from 
its mass m*, radius i?*, and adiabatic index 7 by 
solving the Lane-Emden equation ([7|. 

(3) Use p{r) to determine a distribution of positions X* 
for the fluid elements in Fermi normal coordinates. 



(4) Use Eqs. (13), (14), (h5|, and (B3) to determine 



the positions of the fluid elements in Boyer- 
Lindquist coordinates and their constants of mo- 
tion E + AE, Lz + ALz, and Q + AQ. 

(5) Determine from their positions and constants of 
motion whether each fluid element is on a geodesic 
that plunges directly into the event horizon, es- 
capes to inflnity, or becomes bound to the SBH. 

(6) If a fluid element is o n a b ound orbit, integrate the 
equations of motion (AlO) out to apocenter begin- 
ning from t = 0; the flnal time tf will be half of the 
radial period. 



(7) Create a histogram of the resulting distribution of 
orbital periods to determine the fallback accretion 
rate dm/dt, similar to that depicted in Fig. [2] in the 
Newtonian limit. 



IV. RESULTS 

We begin by presenting our predi ctions for the distri- 
bution d^m/ dEdLzdQ in Sec. IV A this is intermediate 



step (4) in the prescription given above. While this pro- 
vides some insight into how SBH spin affects the tidal 
debris, those interested only in the correspondin g fall- 
back accretion rates dm/dt can skip ahead to Sec. IV B 



Distribution of orbital constants 



Eqs. (13) - (16) reveal that AE, AL^, and AQ for 



each fluid element are inner products of the position X* 
of that fluid element in Fermi normal coordinates with 
covariant derivatives in the direction of the spatial basis 
vectors Af . In the Newtonian limit, these quantities are 
inner products of r — with flat-space gradients. We 
have already seen this in Eq. ([5|, the Newtonian limit of 
Eq. ([l3|. This implies that the projected distributions 
dm/dE^ dm/dLz, and dm/dQ will have the exact same 
shape as the Newtonian distribution shown in Fig.[T] The 
width of the relativistic distribution dm/dE is deflned as 



(17) 



the widths cfl^-, and cfq are found by replacing E with 
Lz or Q in Eq. (17). These widths depend on the SBH 



mass M, dimensionless spin a/M, and the orbit of the 
tidally disrupted star. We explore this dependence in this 
section. 

In the Newtonian limit, (Je ^ E^id given by Eq. ([T]). 
We did not calculate the corresponding widths L^^tid and 
Qtid of dm/dLz and dm/dQ in Sec.jllj as these quantities 
did not enter into the Newtonian orbital period given by 
Eq. ([2|. As the radial period for bound orbits in the Kerr 
spacetime does depend on Lz and Q, we now calculate 
these widths. In the Newtonian limit, 

ALat = (r - Tp) X V 

implying that 
ALz = ALat • z 
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{vp cos L — Lat cos Op) ■ (r — r^) (19) 



and 



AQ 



2 ALat • [(Lat • x)x + (Lat • y)y] 
4GM(rpSin^^ + LArcos6>pCOS^) • (r - r^) .(20) 



SBH 



FIG. 3: The SBH is located at the origin with its spin S point- 
ing along the z axis. The inclination l is the angle between S 
and the star's orbital angular momentum Lat, and Op is the 
angle between S and the position of the star at pericenter. 



where Lat is the orbital angular momentum and the an- 



gles Op and i are defined in Fig. 3 Eqs. \m and m\ 
allow us to define the widths of the Lz and Q distribu- 
tions in the Newtonian limit: 



f2GM\ 



1/2 



^2;,tid 



(cos^^ + cos^(9p)^/^i?* ,(21a) 



Qtid = 4GM(sin'^^ + cos^(9pCos^0^^^^* • (21b) 

Examining the ratio crx/-^tid of the widths of the dis- 
tribution dm/dX of a quantity X in the relativistic and 
Newtonian limits will allow us to assess the accuracy of 
the Newtonian calculation. 



FIG. 4: The widths cfe and (tl^ of the projected distributions 
dm/dE and dm/dLz as functions of the pericenter radius 
Tp/M. These are normalized by the widths ^tid and Lz, tid 
of the corresponding Newtonian distributions. The solid red 
curves are for nonspinning black holes, while the dashed or- 
ange, dot-dashed green, and dotted blue curves are for dimen- 
sionless spins a/M — 0.5, 0.9, and 0.99 respectively. For each 
value of the spin, the curve extending to smaller values of Vp 
corresponds to prograde = 0°) orbits while the other curve 
corresponds to retrograde [l — 180°) orbits. 



We show these ratios for several additional values of 
the dimensionless spin a/M m. Fig. [4j At fixed values 
of Tp, (Je decreases (increases) with a/M on prograde 
(retrograde) orbits. However, at the radius rmb of the 
marginally bound circular orbit (minimum value of 
for each curve), cr£;/£^tid = independent of SBH spin. 
By contrast, the ratio (JL^/Lz^tid at rmb depends sharply 
on a/M. Although CF^^/Lz^tid = 1 all the way down to 
^mb = 4M for nonspinning SBHs, it increases to 2 (de- 
creases to ~ 0.929) as a/M ^ 1 for prograde (retrograde) 
orbits. 



1. Equatorial orbits 

Equatorial orbits {Op = 90°) can be either prograde 
{l = 0°) or retrograde {i = 180°). These orbits have 
Q = 0, and Eq. ( 21b| ) shows that the tidal debris will 
also have Q = in the Newtonian limit. This holds for 
the relativistic calculation as well. If we define 



fr = l 



2GM 



(22) 



2. Inclined orbits 

Inclined orbits have 0° < ^ < 180°; for orbits with 
< it is convenient to define a supplementary incli- 
nation 



(24) 



then for nonspinning SBHs, 

,-1/2 

TTi Jr ^ 

^tid 



COS^ fr^ COS^ Op 



.1/2 



COS^ L + COS^ Orj 



(23a) 
(23b) 



such that l' = 180° — i. If w e define a complementary 
angle Op = 90° — 6>p, Eq. ( AlOc) implies that for parabolic 
{E = 1) orbits 



i < 



-l' < O'p < ^i' 



^ < 90° , 
^ > 90° . 



(25a) 
(25b) 
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FIG. 5: The ratio as /End for orbits with rp = 4M, a/M = 
0.99, and l < 90° as a function of inchnation l and latitude 
Op = 90° — Op. The yellow region above the diagonal is un- 
physical as it corresponds to a portion of parameter space 
forbidden by the inequality (25a). 



The relativistic corrections crx/-^tid to the width of a 
distribution dm/dX will depend on both i and Op for 
inclined orbits, in addition to r^/M and a/M. 

We show (JE/Efid for inclined orbits with Vp = 4M, 
a/M = 0.99, and i < 90° in Fig. [s] The origin of this 
figure corresponds to the same orbit as that shown by 
the point at = 4M on the lower dotted blue curve in 
the upper panel of Fig. [4j We see that the relativistic 
correction CE/Etid is largest when tidal disruption oc- 
curs in the equatorial plane {Op = 0°) and increases with 
inclination t. Although not apparent from this figure, for 
higher spins and smaller radii the ratio (lE/Etid can ex- 
ceed 1.5 on highly inclined orbits, greater than the value 
~ ^/2 found on the marginally bound orbits in Fig. [4] 

The ratio (JE/Etid foi" orbits with Vp = 6M, a/M = 
0.99, and i > 90° is shown in Fig. [6j We must use a 
larger value of than in Fig. [5j as the marginally bound 
orbits with Lz < have larger radii. The origin of this 
figure corresponds to the same orbit as that shown by the 
point at rp = 6M on the upper dotted blue curve in the 
upper panel of Fig. [4j As was the case for l < 90°, the 
relativistic correction is greatest in the equatorial plane 
and increases with inclination l. Figs. [Sland [6] are both 
consistent with our impression from FigT^that the value 
of (JE/Etid largely depends on the difference between rp 
and the radius of the marginally bound orbit for that 
spin and inclination. 

While £^tid is independent of orbital inclination, both 
the numerator and denominator of t he r atio (JL^/Lz^tid 
depend on i and Op. Eq. (21a) and (23b) show that al- 
though both (JL.and Lz^nd ^ as ^ ^ 90°, 6>p ^ 90° 
for a/M = 0, their ratio is undefined in this limit. For 
a/M > 0, the location at which ctl^ =0 shifts to 



FIG. 6: The ratio as/Eud for orbits with rp = 6M, a/M = 
0.99, and l > 90° as a function of supplementary inclination 
l' = 180° — L and latitude Op. The portion of the plot above 
the diagonal is unphysical. 
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FIG. 7: The ratio aLjL^^tid for orbits with rp = 4M, a/M = 
0.99, and l < 85° as a function of inclination l and latitude 
Op. The region of the plot above the diagonal is unphysical. 



^o(a/M) < 90°, O p 9 0° while the Newtonian denomina- 
tor Lz^tid of Eq. (21a) remains unchanged. This implies 
that the ratio (Jl^/ Lz^tid will have a pole at ^ = 90°, 



^p — 0°, as can be seen in Fig. where we show this 
ratio for orbits with rp = AM, a/M = 0.99, and ^ < 85°. 
The figure is cut off above i = 85° to avoid this pole. 
The presence of this pole is also apparent in the contour 
plot of (JL^/Lz^tid for orbits with = 6M, a/M = 0.99, 
and < 85° shown in Fig. [sj We also see that cfl^ ^ ^ 



as one approaches 



80°, 0' = 0° 



Finally we come to the width of distribution of the 
Carter constant Q, which for nonspinning SBHs is given 
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FIG. 8: The ratio ctl, /Lz,tid for orbits with = 6M, a/M = 
0.99, and l > 95° as a function of supplemental inclination 
l' = 180° — L and latitude Op. The region of the plot above 
the diagonal is unphysical. 
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FIG. 10: The ratio aq/Qud for orbits with rp = 6M, a/M = 
0.99, and l > 90° as a function of supplemental inclination 
l' = 180° — L and latitude Op. The region of the plot above 
the diagonal is unphysical. 

(Jq no longer separates into two distinct terms, only one 
of which is dependent on Op. 



B. Fallback accretion rate 

We now present the fallback accretion rate dm/dt for 
stellar tidal disruptions with different values of the SBH 
mass M, SBH spin a/M, pericenter distance r^, orbital 
inclination ^, and polar angle Op at pericenter. We fix 
= Mq, i?* — i?0, and 7 = 5/3 throughout this sec- 
tion, as our focus is on relativistic effects that are inde- 
pendent of these parameters. We measure time in units 
of ttid given by Eq. (|3|, and Tp in units of M (where 
G = c=l). 



FIG. 9: The ratio crg/Qtid for orbits with rp = 4M, a/M = 
0.99, and l < 90° as a function of inclination l and latitude 
Op. The region of the plot above the diagonal is unphysical. 
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We show this ratio for orbits with Tp = 4M, a/M = 0.99, 
and L < 90° in Fig. |9] This contour plot differs qualita- 
tively from the nonspinning result (26) in that crg/Qtid 
is maximized at ^ = 90°, 0^ = 0° instead of near ^ = 0°, 
Op = 0°. We show (Jq/Qtid for orbits with r„ = 6M, 



a/M = 0.99, and i > 90° in Fig. 10 This figure more 



closely resembles the nonspinning result (26), although 



1. Equatorial orbits 

In Fig. [IT] we show the fallback accretion rate dm/dt 
for TDEs of stars on equatorial {Op = 90°), parabolic 
{E 0) orbits with Vp 6M by SBHs with M = lO^M©. 
The different curves correspond to different values of the 
SBH spin and either prograde or retrograde orbits as de- 
scribed in the caption. We see that relativistic correc- 
tions can substantially affect the accretion rate for stars 
with such small pericenters. The peak accretion rate 
dmpeg,]^/dt can be almost twice as large as the Newto- 
nian prediction, and can occur at a time tpeak only half 
that of the Newtonian prediction. SBH spin magnitude 
and direction also strongly affect the accretion rate; large 
prograde = 0°) spins reduce dm^Q^\^/dt and delay tpeak, 
while large retrograde spins {i = 180°) have the opposite 
effect. The jaggedness in the curves results from the fi- 
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FIG. 11: The dimensionless fallback accretion rate 
{t tid / m^) dm / dt as a function of time t/ttid following the tidal 
disruption of a stars on equatorial orbits with Vp = 6M by 
SBHs with M = 10^ M0. The smooth, solid black curve is the 
Newtonian prediction shown in Fig. |2] The solid red curve is 
the relativistic prediction for a non-spinning SBH. The dashed 
orange (dotted blue) curves correspond to spin magnitudes 
a/M = 0.5 (0.99) respectively. The curves of each type with 
lower (higher) peak accretion rates correspond to orbital in- 
clinations L = 0° (180°). 



FIG. 12: The dimensionless fallback accretion rate 
{ttid/TTi^)dm/dt as a function of time t/tud following the tidal 
disruption of a stars on equatorial orbits with Vp = 6M by 
SBHs with M = 10^ M©. The smooth, solid black curve is the 
Newtonian prediction. The solid red curve is the relativis- 
tic prediction for a non-spinning SBH. The dashed orange, 
dot-dashed green, and dotted blue curves correspond to spin 
magnitudes a/M of 0.5, 0.9, and 0.99 respectively. The curves 
of each type with higher (lower) accretion rates at late times 
correspond to orbital inclinations ^ = 0° (180°). 



nite number of particles {N = 10^) and finite proper time 
step (Ar = lO^M). Each curve took approximately 15 
minutes to prepare on a single-processor 2.2 GHz lap- 
top; better optimization and additional computational 
resources could allow a more thorough exploration of pa- 
rameter space with greater accuracy. 

In Fig. [T2j we again show the fallback accretion rate 
for equatorial orbits with = 6M, but have reduced the 
SBH mass from lO^M© to lO^M©. By measuring time in 
units of ttid, we have eliminated the the dependence of the 
Newtonian prediction on M and r^. The magnitude of 
the relativistic effects is similar to that in Fig. [TT] for the 
prograde orbits, but the accretion rate has been sharply 
cut off at early times for the retrograde orbits. This cutoff 
occurs because the mostly tightly bound particles (those 
which would have the shortest orbital periods) have lost 
so much energy E and angular momentum Lz that they 
are directly swallowed by the SBH event horizon. This 
effect is very small for a/M = 0.5, where only 0.021% 
of the star's initial mass is directly swallowed (a slight 
reduction in the retrograde dashed orange curve at the 
earliest times is just barely noticeable). The cutoff is 
much more pronounced for a/M = 0.9 and 0.99 shown by 
the dot-dashed green and dotted blue curves respectively. 
In these cases 21.4% and 34.2% of the initial stellar mass 
is directly captured. This directly captured material will 
have little opportunity to radiate, and should therefore 
not contribute to observed TDE light curves. 

In our dimensionless units, the Eddington accretion 



rate is 



ttid dm 



Edd 



dt 



An'^mp 
r]ml(jTC' 



-3/2 



M 



(gm) (lO^Mcr 



7/2 



(27) 



where is the proton mass, r] is the radiative efficiency, 
and (Jt is the Thomson cross section. While the peak 



accretion rates for 10 M© SBHs as shown in Fig. 11 
are comfortably sub-Eddington, those for 10^ Mq SBHs 
as shown in Fig. [12] are highly super-Eddington. It is 
therefore unlikely that the observed TDE luminosity for 
M = lO^M© will track the fallback accretion rate at 
early times which are most sensitive to SBH spin. More 
massive SBHs have peak accretion rates below the Ed- 
dington limit, but are only capable of tidally disrupting 
stars when 



M ^ M 







-2/3 



Mo 



-1/3 



(28) 

SBHs as massive as 10^ M© will directly capture many 
stars without tidally disrupting them, although relativis- 
tic effects allow TDEs for M < lO^M© despite this New- 
tonian estimate ^38n40 . These considerations suggest 
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FIG. 13: The dimensionless fallback accretion rate 
{ttid/rn^)dm/dt as a function of time t/tud following the tidal 
disruption of a stars on equatorial orbits with = 12M by 
SBHs with M = 10^ Mq. The smooth, solid black curve IS 
the Newtonian prediction. The solid red curve is the rela- 
tivistic prediction for a no n- spinning SBH. The dashed or- 
ange (dotted blue) curves correspond to spin magnitudes 
a/M of 0.5 (0.99). The curves of each type with lower 
(higher) peak accretion rates correspond to orbital inclina- 
tions ^ = 0° (180°). 



FIG. 14: The dimensionless fallback accretion rate 
{ttid/rn^)dm/dt as a function of time t/ttid following the tidal 
disruption of a stars on inclined orbits with Vp = 6M by SBHs 
with M = 10^ Mo and a/M = 0.99. The smooth, sohd black 
curve is the Newtonian prediction. The remaining curves cor- 
respond to relativistic predictions with polar angle Op = 90° 
but different values of the orbital inclination t. The peak ac- 
cretion rate drripeak/dt increases with l; with increasing values 
of drupeak/dt, the dotted blue, dot-dashed green, dashed or- 
ange, solid red, dotted purple, and dot-dashed cyan curves 
correspond to ^ = 0°, 45°, 75°, 105°, 135°, and 180°. 



that M :^ lO^M© may be a sweet spot for attempts to 
constrain SBH spins with TDE light curves, which is for- 
tuitous since SBHs of this mass are fairly ubiquitous at 
the modest redshifts where TDEs are readily observable. 

We show the fallback accretion rates for TDEs by a 
SBH with M = 10^ Mq in Fig. [l3) The pericenter for 
the TDEs shown in this figure has been increased to 
Tp = 12M to show how relativistic effects fall off with dis- 
tance. Comparing this figure to Fig. [TTJ we see that rel- 
ativistic corrections increase dm^QQ\^/dt and reduce tpeak 
by at most ~ 25%, and SBH spin induces similar correc- 
tions to these quantities at the 25% level. 



2. Inclined orbits 

Tidally disrupted stars are scattered onto their initial 
orbits at large radii where the influence of the SBH spin is 
negligible. One should therefore expect a flat distribution 
in cos^ in the range — 1 < cos^ < 1. We now consider 
how the fallback accretion rate depends on l and the polar 
angle Op at pericenter for inclined orbits. We show dm/dt 
for fixed spin magnitude a/M = 0.99 and O^j = 90° but 



different orbital inclinations l in Fig. 14 The Newtonian 
prediction in this figure is identical to that in Fig. [TTJ and 
the dotted blue (dot-dashed cyan) curves corresponding 
to ^ = 0° (180°) are identical to the dotted blue curves 
m Fig. [11] The remaining curves correspond to inclined 



orbits that were discussed in detail in Sec. IIV A 21 In 
terms of the relativistic correction crE/End to the width 
of the energy distribution, increasing t from 0° to 180° 
while keeping Op = 90° fixed corresponds to moving from 
left to right along the x-axis of Fig. [5j then moving from 
right to left along the x-axis of Fig.jo] The ratio (TE/Etid 
increases with l along this path, qualitatively implying 
that tidal debris falls deeper into the SBH's potential 
well with shorter radial periods t/M. Although in the 
relativistic limit the radial period depends on Lz and 
Q in addition to consideration of E alone seems to 
illustrate the qualitative dependence on t. 

We show the dependence of the fallback accretion rate 
on the polar angle at pericenter Op in Fig. 15 As in 



the previous figure. Fig. [15] depicts TDEs of stars with 



pericenter Tp 



6M by a SBH with M = 10^ M© and 



a/M = 0.99. The orbital inclination is fixed at ^ = 80°, 
but the latitude O'p = 90° -6>p is varied from 0° to 60°. In 
terms of the ratio (JE/End^ this corresponds to moving 
upwards along the vertical line i = 80° in Fig. [5] The 
ratio (JE/Eti^ decreases as we move along this path, im- 
plying that the tidal debris becomes less tightly bound to 
the SBH with longer radial periods t/M. This qualita- 
tively accounts for the trend seen in Fig. [Tsj that dm/dt 
is shifted to later times as Op decreases (^^ increases). 
Although not shown, this same trend holds for l > 90° 
as one would expect from consideration of Fig. [Oj Our 
results for the peak accretion rate drupeak/dt and the 
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^/^tid 



FIG. 15: The dimensionless fallback accretion rate 
{t tid / m^) dm / dt as a function of time t/ttid following the tidal 
disruption of a stars on inclined orbits with Vp = 6M by SBHs 
with M = lO^Mo and a/M = 0.99. The smooth, sohd black 
curve is the Newtonian prediction. The remaining curves cor- 
respond to relativistic predictions with inclination l = 80° 
but different values of the polar angle at pericenter Op. The 
peak accretion rate drupeak / dt increases with Op ; with increas- 
ing values of drupeak/dt, the dotted blue, dashed orange, and 
solid red, curves correspond to Op = 30°, 60°, and 90°. 



time tpeak at which this peak occurs are summarized in 
Table III 



M 


a/M 


Tp 




Op 


^peak 


ttid ^^peak 
m* dt 




0.0 


6M 


0° 


90° 


1.87 


0.0683 


10^ Mo 


0.5 


6M 


0° 


90° 


2.04 


0.0645 


10^ Mo 


0.5 


6M 


180° 


90° 


1.72 


0.0766 


10^ Mo 


0.99 


6M 


0° 


90° 


2.17 


0.0596 


10^ Mo 


0.99 


6M 


180° 


90° 


1.56 


0.0838 


lO^Mo 


0.0 


6M 


0° 


90° 


1.92 


0.0685 


lO^Mo 


0.5 


6M 


0° 


90° 


2.11 


0.0638 


lO^Mo 


0.5 


6M 


180° 


90° 


1.81 


0.0747 


lO^Mo 


0.9 


6M 


0° 


90° 


2.24 


0.0595 


10^ Mo 


0.9 


6M 


180° 


90° 


4.73 


0.0335 


lO^Mo 


0.99 


6M 


0° 


90° 


2.24 


0.0585 


lO^Mo 


0.99 


6M 


180° 


90° 


12.33 


0.0079 


lO^Mo 


0.0 


12M 


0° 


90° 


2.25 


0.0583 


lO^Mo 


0.5 


12M 


0° 


90° 


2.28 


0.0564 


lO^Mo 


0.5 


12M 


180° 


90° 


2.17 


0.0599 


lO^Mo 


0.99 


12M 


0° 


90° 


2.31 


0.0558 


lO^Mo 


0.99 


12M 


180° 


90° 


2.05 


0.0617 


10^ Mo 


0.99 


6M 


45° 


90° 


2.09 


0.0636 


10^ Mo 


0.99 


6M 


75° 


90° 


1.89 


0.0665 


10^ Mo 


0.99 


6M 


105° 


90° 


1.78 


0.0735 


10^ Mo 


0.99 


6M 


135° 


90° 


1.61 


0.0781 


10^ Mo 


0.99 


6M 


80° 


30° 


2.15 


0.0606 


10^ Mo 


0.99 


6M 


80° 


60° 


2.02 


0.0651 


10^ Mo 


0.99 


6M 


80° 


90° 


1.89 


0.0680 



TABLE I: Peak accretion rate drripeak/dt and the time tpeak 
at which this peak occurs for different values of the SBH mass 
M, spin a/M, pericenter r^, orbital inclination and polar 
angle Op at pericenter. 



V. DISCUSSION 

Our primary goals in this paper have been to estab- 
lish a theoretical framework for investigating stellar tidal 
disruption deep in the relativistic regime and to conduct 
a preliminary survey of the many dimensional parameter 
space associated with this problem. Our starting point 
for this framework was the Newtonian model proposed in 
LKP09, in which a star's orbit and density profile were 
used to determine the distribution of orbital energy of 
the tidal debris and hence the fallback accretion rate. 
The gravitational potential of a Newtonian point mass is 
spherically symmetric, implying that for a given stellar 
profile the only physically significant parameters are the 
SBH mass M and the pericenter Tp. The Kerr space- 
time however is only axisymmetric, increasing the num- 
ber of required orbital parameters to include the SBH 
spin a/M, the orbital inclination ^, and the polar angle 
at pericenter Op. The radial period of the tidal debris 
now depends on the angular momentum and Carter 
constant Q in addition to the energy E. Fermi normal co- 
ordinates allow us to determine the distribution of these 
orbital constants for a given stellar profile and orbit, and 



the Killing vectors and tensor of the Kerr spacetime pro- 
vide a series of first-order differential equations for the 
orbital motion of the tidal debris given the orbital con- 
stants. Combining these ingredients according to the pre- 
scription given at the end of Sec.jTTTj we can calculate the 
fallback accretion rate dm/dt as a function of time t. 

If dm/dt is greater than the Eddington accretion rate, 
as will be the case for M < l{f Mq and t ~ tpeak, the 
emission may be dominated by powerful super-Eddington 
outflows whose luminosity will not be proportional to 
dm/dt |53jj. At later times and for larger SBH masses, 
the luminosity will be dominated by an accretion disk 
about the SBH, whose emission may be partially repro- 
cessed by unbound tidal debris. The bolometric lumi- 
nosity of this accretion disk should be proportional to 
dm/dt, although the finite extent of the disk implies that 
the emission spans a fairly narrow range of temperatures 
j45] . Monochromatic light curves may therefore not be 
proportional to dm/dt; optical emission will lie on the 
Ray leigh- Jeans tail of the distribution and thus scale as 
{dm/dty^^ [ill [53l. Although the connection between 
observed monochromatic light curves and the fallback ac- 
cretion rate is highly nontrivial, this accretion rate con- 
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stitutes the primary input for more sophisticated calcu- 
lations of TDE emission. 

Our preliminary survey of relativistic tidal disruption 
suggests that the spin magnitude a/M, inclination ^, and 
the polar angle at pericenter 6p can all affect the fall- 
back accretion rate provided is sufficiently small. The 
newly discovered optical transient PSl-lOjh [13 may cor- 
respond to a TDE with = 12M [54 where relativistic 
effects can be significant as seen in Fig. 13 The light 



curve associated with this event shows systematic differ- 
ences with the LKP09 model during the early rise that 
may be sensitive to relativistic corrections; a detailed 
comparison between our model and this light curve is 
an important subject of future work. 

We find that as approaches the radius of the 
marginally bound circular orbit (below which the star 
will be directly captured by the SBH), relativistic cor- 
rections can increase the peak accretion rate dmpeg,\^/dt 
and reduce the time tpeak at which this peak accretion 
occurs by a factor of two. For a fixed value of r^, rel- 
ativistic effects are largest for large SBH spins that are 
anti-aligned with the star's initial orbital angular mo- 
mentum. However, as explored in detail in Sec. |IVB[ 
there is a great deal of degeneracy between a/M, ^, 
and 6p. This degeneracy is not terribly surprising given 
the crudeness of our model, as each of these parame- 
ters can only increase or decrease the distribution of fall- 
back times. Multi-frequency observations may be able to 



partially break this degeneracy. Comparison of Figs. 11 



and [14] reveals a degeneracy between the fallback accre- 
tion rates for TDEs of stars initially on inclined orbits 
and those on equatorial orbits of SBHs with lower spins. 
Although the fallback accretion rates may be the same, 
the disks that subsequently form about the more highly 
spinning SBHs should have smaller radii and thus higher 
temperatures. Emission from these disks should there- 
fore be harder than that from the disks of SBHs with 
lower spins. A more quantitative study of the degener- 
acy between relativistic parameters is another subject of 
potential future work. 

Our model of tidal disruption, like that of LKP09, is 
an example of what Guillochon and Ramirez-Ruiz 
describe as a "freezing" model in that the distribution of 
the orbital constants of the tidal debris is "frozen in" at 
pericenter. This failure to account for the redistribution 
of orbital constants during the small but nonzero time 
over which tidal disruption occurs can lead to several 
discrepancies in the fallback accretion rate as discussed 
extensively in Guillochon and Ramirez-Ruiz [55 . For 
fully-disruptive encounters like those considered in this 
paper, hydrodynamical simulations performed by Guillo- 
chon and Ramirez-Ruiz suggest that dm/dt is insensitive 
to Tp when it is much below the tidal radius rtid, in con- 
trast to the predictions of freezing models. The star's 
self-gravity is insufficient to hold it together at pericen- 
ter for rp <C rtid, so it is inappropriate to apply freez- 
ing models to the unperturbed star at Tp in this regime. 
Future work could explore applying relativistic freezing 



models to the star at rtid instead of Tp as suggested by 
Guillochon and Ramirez-Ruiz; we would expect relativis- 
tic corrections to dm/dt to be much smaller in this case. 

Another problem of freezing models is that they fail to 
account for pressure gradients that develop when the star 
is tidally compressed [48l [56] . These pressure gradients 
can redistribute material to more tightly bound orbits, 
leading to an increased feeding rate at early times. Our 
study suggests that this increase in the early rise could 
be misinterpreted as a relativistic effect if Vp is small. We 
hope to go beyond the limitations of freezing models in 
future work on relativistic tidal disruption. One subject 
deserving special scrutiny is the conditions under which 
the tidal debris viscously evolves into an accretion disk 
about the SBH. Lense-Thirring precession [35 should de- 
lay this process for stars on inclined orbits, helping to 
further break some of the degeneracy between spin mag- 
nitude a/M and inclination l. 

Despite this need for additional work, the future 
prospects for using TDEs to measure SBH spins seem 
highly promising. The number of observed TDE candi- 
dates has increased dramatically in the past few years, 
and future transient surveys should provide hundreds if 
not thousands of additional candidates [M]. A portion 
of these candidates, like PSl-lOjh [13], should have suffi- 
ciently small Tp such that relativistic effects will be signif- 
icant. We look forward to challenging our model and its 
future refinements with this rich observational bounty. 
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Appendix A: Kerr metric 

In this appendix, we provide explicit formulae for 
quantities relevant to describing tidal disruption in the 
Kerr spacetime [57 . In Boyer-Lindquist coordinates [58^ 
and units where G = c = 1, the Kerr metric takes the 
form 



ds^ = ga(3dx^dx^ 



2Mr\ , ^ 4Marsin^6> , , , 
1 — dt^ dtdcj) 

—dr^ + S ^6*2 + ^ sin^ 6» #2 



(Al) 



where M is the mass of the black hole, a/M is its dimen- 
sionless spin, and 



S = r^ + cos^ , 

A r^ + - 2Mr , 

A = {r^ ^ a^f - Aa^ sinh 



(A2a) 
(A2b) 
(A2c) 
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The inverse of this metric is given by 
/ 
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(A3) 
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Massive particles move on timelike geodesies of the Kerr 
metric with 4-velocity 



{i f 



(A4) 



where overdots denote derivatives with respect to proper 
time r. The Kerr metric possesses timelike and azimuthal 
Killing vector fields, which we will denote by {d/dt)'^ and 
{d/dcj))^ following the notation of Wald [52 . This implies 
the existence of a conserved specific energy and angular 
momentum 



E 



dt 



/3 



d 



/3 



(A5a) 



(A5b) 



The Kerr metric also possesses a Killing tensor [521 159] 

Ko,(3 = S(/a^/3 + nJis) + r'^ga/3 , (A6) 



where 

r = 



2S 



a 
A 

a 

2S 



d 



f Q ^ 


a 




) 




( d 


-( 

2S ^ 


\ dr 



(A7a) 

a 

(A7b) 



are the repeated principal null vectors found in a Petrov 
classification of the Weyl tensor [S2l EQI [611 . This implies 
the existence of a third constant of motion 



K 



(A8) 

The Carter constant Q is defined in this way so that in the 
Newtonian limit (r ^ oo, E ^ 1), Q ^ L^. ^ L'^^ where 
Lx and Ly are the x and y components of the orbital 
angular momentum. We can also define the inclination 

cos i , ( A9) 

which in the Newtonian Hmit converges to the angle be- 
tween t he S BH sp in a nd angular momentum Lyv- 

Eqs. ( |A5| ) and (A8), along with the normalization of 
the 4-velocity gapu°'u^ = — 1, can be rearranged to pro- 
vide equations of motion on Kerr geodesies [62] 

AE-2MarL^ 

X (AlOa) 



aL,f 
- aEf 



■Q] 



Q 

Lz csc^ 6 



L\ cot^ i9 - a^(l -£;^)cos^ 



2MarE - a^L^ 



(AlOb) 
(AlOc) 

(AlOd) 



Evaluating the covariant derivatives of tensors in 
Boyer-Lindquist coordinates requires the use of Christof- 
fel symbols F^^ [52] 



a-'- ci-'-ci — ^a-'- ci---( 



•d---bk 



.^rtr"--"^ c,...d...c, .(All) 



The Christoffel symbols are symmetric under exchange of 
their lower two indices. In Boyer-Lindquist coordinates, 
the non-vanishing Christoffel symbols are 

-(2r^ - S) 



^ tt 



^ tt 



J- ee 
J- ee 

j.e 
^t 



te 



^ tr 

J.t 
J- te 

J- tcj) 

j^e 
T^e 

J- 



r* 



rcf) 



rcf 

^e^ 
p0 



S3 

2Ma^rsin6>cos6> 

S3 

r r - M 



a" sm U cos ( 



EA 



rA 

sin ^ cos t 
A sin^ (9 
sin 6> cos 



Ma^ sin^ 

S2 



(2r' - S) 



(A12a) 
(A12b) 
(A12c) 
(A12d) 
(A12e) 
(A12f) 
(A12g) 



S3 



[(r^ + a^)A - SAa^ sin^ i9] (A12h) 



S^A 



(2r" - S) 



Ma 



(2r' - S) 



2Ma^r sin l9 cos (9 

2Mar cos 6> 

sin^ 

MaA sin^ 6> , ^ 
(2^'-^) 

2Mar{r'^ + a^) sin 6> cos 6> 

S3 

sin 6> cos 



r 
S 



Ma sin^ 



SA 



2r^ 



(r^ + a^) + r' 



^ a^ 



a^ sin^ < 



M 



2Mr^ 



S SA 
2Ma3rsin^6>cos6> 

COS l9 / 2Ma^rsin^l9 
1+ 



(A12i) 

(A12j) 

(A12k) 

(A121) 

(A12m) 

(A12n) 

(A12o) 
(A12p) 

(A12q) 

(A12r) 

(A12s) 
(A12t) 
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Appendix B: Fermi normal coordinates 



where 



Although the Boyer-Lindquist coordinates utilized in 
the previous appendix are an excellent choice for describ- 
ing motion along geodesies of the Kerr spacetime, Fermi 
normal coordinates [63] provide a better description of 
the local neighborhood of a particle (such as the center 
of mass of a tidally disrupting star) moving along such 
a geodesic. These coordinates are defined by first choos- 
ing a reference point Pq on a central geodesic G as the 
origin. One then chooses an orthonormal tetrad of 4- 
vectors in the tangent space at Pq? where Aq is the 
tangent vector to G and A^ {i = 1, 2, 3) are three space- 
like vectors. The point (r, X*) (i — 1,2,3) in Fermi nor- 
mal coordinates is reached by traveling a proper time r 
along G to a point h(T\ then traveling a proper distance 
s = [Xl^(-^*)^]^^^ along the orthogonal geodesic passing 
through h{r) whose tangent vector is X^A^. 

Marck [64 constructed an explicit tetrad of 4- vectors 
A^ that are parallel propagated along timelike geodesies 
of the Kerr spacetime. He provided components of these 
4- vectors in the canonical symmetric orthornormal tetrad 
first introduced by Carter [65], but to avoid introducing 
yet a third set of coordinates we provide these vectors in 
Boyer-Lindquist coordinates. Aq is simply the 4- vel ocity 
of the center of mass, given explicitly by Eqs. (A4) and 



(AlO). The Boyer-Lindquist components of the spacelike 



4- vectors A? are 



A? 



ar 



a{r'^ + a?)rf 



A 

„2 I „2\ 



(5a sin 9 cos ( 



a I arr /3cos( 



4 



AS 



ifl/2 



smt 



(r^ + a^) cos^f 



A 

„2 I „2\ 



rsmt 



acos0 ; .,N , 
[E{r'+a')-aL,] 



A = 



r 



Sifl/2 



aE sin 9 



1 



a? cos Or rO 



/fi/2 



smf 



(Bla) 
(Bib) 
(Blc) 

(Bid) 

(Ble) 
(Blf) 
(Big) 

(Blh) 



a 



K-a" cos^^^'/' 



1 

a 



(B2a) 
(B2b) 



The basis 4-vectors Ai and A3 given above are equal to 
Ai and A3 given in [64l, and are not parallel propagated 
along the central geodesic G. However, since we only 
make use of Fermi normal coordinates at a single point 
along G (the pericenter of the tidally disrupted star's or- 
bit), we are justified in using these vectors in our tetrad. 
If the origin Pq (the center of mass of a star) is located 
at Xq in Boyer-Lindquist coordinates, an event located 
at (0,X*) in Fermi normal coordinates will be located at 



(B3) 



AS 



= a- 



xt = 



EA 



-[£(r2 + a2)-aL,] 



-p-{aEsm^9-L,) 
ar 
P9 

^[E{r' + a')-aL.]-^{aE- 



(Bli) 

(Blj) 
(Blk) 

(Bll) in Boyer-Lindquist coordinates. 
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